arXiv: 1506.0345 lv2 [astro-ph.GA] 21 Aug 2015 


Mon. Not. R. Astron. Soc. 000. [THTSl 0 Printed 24 August 2015 (MN IATfrjX style file v2.2) 


Carbon-enhanced metal-poor stars in dwarf galaxies 


Stefania Salvadori 1 ’*, Asa Skiiladottir 1 & Eline Tolstoy 1 

1 Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 97f 7 AD Groningen, The Netherlands 
-kVENI Fellow 


ABSTRACT 

We investigate the frequency and origin of carbon-enhanced metal-poor (CEMP) stars 
in Local Group dwarf galaxies by means of a statistical, data-calibrated cosmological 
model for the hierarchical build-up of the Milky Way and its dwarf satellites. The 
model self-consistently explains the variation with dwarf galaxy luminosity of the 
observed: i) frequency and [Fe/H] range of CEMP stars; ii) metallicity distribution 
functions; iii) star formation histories. We show that if primordial faint supernovae 
dominated the early metal enrichment, then CEMP-no stars enriched by the first 
stellar generations should be present in all dwarf galaxies, with similar number of 
stars and CEMP fractions at [Fe/H]< —4. We demonstrate that the probability to 
observe a star that is carbon-enhanced within a given [Fe/H] range strongly depends 
on the luminosity of the dwarf galaxy and, on average, it is an order of magnitude lower 
in “classical” Sculptor-like dSph galaxies (P ^ 0.02) than in the least luminous ultra¬ 
faint dwarfs (P « 0.1). In addition, we explain why it may be easier to find CEMP-no 
stars at [Fe/H]ss —2 in classical dSph galaxies than in ultra-faint dwarfs. These are 
consequences of the dramatic variation in the fraction of stars at [Fe/H]< —3 with 
galaxy luminosity: ^ 40% for galaxies with L < 10 5 Lq, and < 0.2% for L > 10 7 Lq. 
We present model predictions for the low Fe-tail and CEMP fraction of stars in dwarf 
galaxies, with particular emphasis on the Sculptor dSph, that can be used to shed 
light on the properties of the first stars. 

Key words: galaxies: high-redshift, dwarf, Local Group; stars: abundances; cosmol¬ 
ogy: theory. 


1 INTRODUCTION 

Understanding the nature of the first stars is one of the 
key questions of modern Cosmology. Despite extensive ob¬ 
servational searches, truly pristine stars have so far escaped 
detection. Currently, the most metal-deficien t star kn own 
has a total metallicity of Z « 6.9 x 10 ' dCaffau et al.l 
l201ll l. The persistent lack of more metal-poor stars seems 
to confirm the idea tha t primordial stars were all more mas- 
sive t han ~ 1M& (e.g. McKee fe Tanll2008l : iHosokawa et al.l 
l201ll : iHirano et al.ll2014l ), and that th eir formation was pos¬ 
sibly quenched at ea rly times (e.g. ISalvadori et al.l 120071 : 
IPallottini et al.l I2014T ). However, the chemical signatures 
of these extinct stellar generations could be retained in 
the photospheres of ancient (> 12 Gyr) low-mass second- 
generation stars, which formed in pre-enriched environ¬ 
ments, Z > Z cr = 10 ~ 5±1 Zq, where metals and dust grains 
dispersed by the fir st stars enab led efficie nt gas cooling and 
fragmentation (e.g. ISchneider et alj|2003 l. These stellar fos¬ 
sils should be observable in the oldest stellar components of 
our Galaxy and in its ancient dwarf satellites. 

High- and medium-resolution spectroscopic studies 
of Galactic halo stars have revealed the existence of a 


popu l ation of carbon-rich stars (e.g. [Beers_&_Chnstlie3 
20051: lAoki et al.ll2007l : lYong et aHl2013l : iNorris el al.ll20lT 


Lee et alJ I20TT . These objects are u sually defined t o 


have carbon-to-iron ratio [C/Fe]> 0.7 (lAoki et alJ 120071 ') . 
and they can be divided into two main populations: 
carbon-rich stars that exhibit an excess in heavy ele¬ 
ments formed by slow (or rapid) neutron capture pro¬ 
cesses, CEMP-s (CEMP-r) stars, and carbon-rich stars 
that do not exhibit such an excess, CEMP-no stars. The 
available data are consistent with the idea that C EMP- 
s stars belong to bina ry systems dLucatello et al.l 120051 : 
IStarkenburg et al.ll2014r i. and have acquired their carbon- 
excess f rom an Asymptotic Giant Branch (AGB) compan¬ 
ion (e .g. lBisterzo et al.l2012l : |placco et ahllSOl .'ll : I Abate et al.l 
l2015ll . Thus, the chemical abundances measured in these 
stars are not representative of the interstellar medium (ISM) 
out of which they formed. On the other hand, CEMP-no 
stars are not preferentially associated to binary systems 
(lHansen et al.l 120131 : Cohen et al.l 120131 : INorris et all 120131 : 
IStarkenburg et ahl l2014l h Hence there is no observational 
evidence supporting the idea of mass transfer as the ori¬ 
gin of their chemical abundances, which was suggested by 
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some authors (e.g. ISuda et all 12004 iKomiva et al.l 120071) . 
Furthermore, both their frequency and carbon-exce ss in¬ 
crease wi t h decreasing [Fe/Hl (e.g . iLucatello et al.l 120051 : 
I Lee et ah ] 120131 : I Norris et al.l 120131 ). and 8 out of the 9 


(Christlieb et al.l 

2002; Frebel et al. 2005; Norris et al. 2007: 

Keller et al. 2014 

: Hansen et al. 20151: Bonifacio et al. 2015 

Allende Prieto et al. 20151: Frebel et al. 20151). These find- 


ings favor the idea that CEMP-no stars are a peculiar stellar 
population and that their chemical abundances likely reflect 
their birth environment. 

The unusual chemical compositions of the most iron- 
poor and carbon-rich stars can be successfully matched 
by models of primordial faint SN that experienced 
mixing and fallback, hence releasing small amounts of 
iron and l arge amounts of carbon and other light ele¬ 


ments (e.g. UmedafcMjomoto 2003; Iwamoto_et_alJ_ 2005J; 


Joggerst et al. 120091 : Marassi et al.l I20l4 : Tominaga et al.l 


2014). A reatively good agreement with observations is 


also obtained by models of zero- (or very low-) metallic- 
ity massive “spinstars”, which experience mixing and mass 
loss because of the i r very high-rotational vel ocities (e.g. 
iMevnet et al.l 1200(1 l20ld : iMaeder et al.l 120151 ). Recently, 
chemical evolution studies have further support the idea of 
the link between primordial faint SN and CEMP-no stars, 
showing that the observed fraction of carbon-enhanced to 
carbon-normal stars at [Fe/H]< —3 is successfully repro- 
duced if faint SN dominated the early metal-e nrichment 
rtde Bennassuti et al.ll20l4 : lCooke fe Madaul20l4 ). Thus, we 
can work under this simple hypothesis to predict the fre¬ 
quency of pristine carbon-enhanced stars in ancient and 
metal-poor dwarf galaxies. 

Carbon-enhanced metal-poor stars have been found 
in a significant fraction in the faintest satellites of the 
Milky Way, the so-called ultra-faint d warf galaxies , wit h 


Frebel et al. 

UOILICO D ^ 1U I^UIIIO Cl CXI. 

2010b!; Lai et al. 201 ll: Gilmore et al. 

2013; 

Frebel et al. 

galaxies to 
minihaloes, 
prior the en 

2014). Different grouns have nronosed these 
be the living relics of the first star-forming 
which formed in the Milky Way environment 
d of reionization (e.e. Salvadori & Ferrara 2009; 

Bovill & Ricotti 2009): Munoz et al. 

2009; Bovill & Ricotti 

20lIF Salvadori & Ferraral 20l"3; 

Salvadori et al. 20l4). 

These theoretical predictions are consistent with recent ob¬ 
servations of star formation histories in ultra-faint dwarf 
galaxies (e.g. Dall’Ora et al. 20121: Okamoto et al. 20121: 


I Brown et al.l 


2014 ). In particular, by interpreting the 


observed Fe-Luminosity relation and metallicity distribution 
function (MDF) of dwarf g alaxies in a cosmological context, 
ISalvadori fc Ferraral (l2009h predicted these ancient ultra¬ 
faint dwarf galaxies to be the best objects to look for the 
chemical imprints of the first stellar generations. This pic¬ 
ture is supported by the recent di scovery of several CEMP- 
no stars at [Fe/H]< —3 in Segue 1 (iFrebel et al.l2014l) . which 
is one of the faintest ultra-faint dwarfs. 

However, carbon-enhanced metal-poor stars seem to 
be rare in the more luminous “classical” dwarf spheroidal 
(dSph) galaxies, L > 10 5 Lq, where measurements are 
available for larger stellar samples. The deficiency of 
CEMP(-no) stars is especially mysterious in the Sculp¬ 
tor dSph galaxy. The observed Color-Magnitude-Diagram 
(CMD) and MDF of Sculptor are consistent with this 


galaxy being dominate d by ancient stars, > 10 Gyr 

old dde Boer et al.l 120121 ). However, no CEMP stars have 
been found among the ten carefully studied stars at 
[Fe/H[< —3 dFrebel et all 2010a: Tafelmever et al. 2010|; 


Starkenburg et all 20131 : Jablonka et al.1 20151 : Simon et al.1 

20151) . Only recently, the first CEMP -no star ha s been dis¬ 


covered at [Fe/H]« —2 dSkuladottir et aTll2015l ). The ob¬ 
served chemical abundance pattern of this star is consistent 
with the idea that the carbon excess is the result of a pris¬ 
tine population of faint SN polluting its birth-environment. 
However, the lack of CEMP-no stars at lower [Fe/H] is not. 
The important question is, are carbon-enhanced metal-poor 
stars missing in classical dSph galaxies, or are they hidden? 

We argue they may be hidden. To sustain this thesis, 
in this paper we use a statistical, data-calibrated cosmolog¬ 
ical model for the formation of the Milky Way and its dwarf 
satellites, to predict the frequency of CEMP (- 110 ) stars in 
Sculptor and other Local Group dwarf galaxies. We present 
a simple, global scenario that self-consistently explains the 
variation with galaxy luminosity of the observed: i) carbon- 
enhanced to carbon-normal star ratios; ii) metallicity distri¬ 
bution functions; and iii) star formation histories. We link 
the properties of Local Group dwarf galaxies with first stel¬ 
lar generations and early galaxy formation processes, and 
present model predictions aimed at identifying the hidden 
CEMP-no stars in the classical dSph galaxy Sculptor. 

For _ consistency with previous work s 

llSalvadori &; Ferraral 120091 . 120121 : ISalvadori et al.l 12014 ) . 
along with Via Lactea II, Aquariu s, and CLUE S sim¬ 
ulations (e.g. Mada.u et al.l [20081 : IWang et al.l 120121 : 
iBem'tez-Llambav et al.l 120151 ). we adopt a Lambda Cold 
Dark Matter (ACDM) cosmology with = 0.24, 

Da = 1 — D m = 0.76, fib = 0.04, and Ho = 73km/s/Mpc. 
Furthermore, in the overall paper we a ssume the solar 
abundance values bv lAsplund et aD J2009I ). 


2 OBSERVATIONS OF CARBON-RICH STARS 

Current observations of CEMP stars, [C/Fe]> 0.7, are shown 
in Fig.l, where we present a sample of stars from the liter¬ 
ature with measured carbon-to-iron ratio, [C/Fe], and iron- 
abundance, [Fe/H], in both the Galactic halo and nearby 
dwarf galaxies. Most of these data are based on one dimen¬ 
sional local thermodynamic equilibrium (1D/LTE) stellar 
model analysis. _ 

Measurements for halo stars are from IPlacco et al.l 
(12014) . who selected among all available data (see references 
therein) a sample of 505 stars with [Fe/H]< —2 and [C/Fe] 
measurements. This sample i nclude s dwarf and giant stars 
(0 < log g < 5). [Placco et al.l (I20l4) corrected [C/Fe] to ac¬ 
count for the depletion of surface carbon abundance, which 
is expected to occur on the upper Red Giant Branch (RGB), 
log g < 2. The correction depends on several observed quan¬ 
tities: logg , [C/Fe], [Fe/H] and also, to a lesser extent, on 
[N/Fe]. Hence, we should keep in mind that the corrected 
values presented in Fig. 1 have intrinsic errors, which are 
estimated by the authors to be always within the 2<r uncer¬ 
tainties of the [C/Fe] measurements, i.e. ±0.3 dex. 

In more distant dwarf galaxies only the brighter RGB 
stars are typically observed. Hence, in our selection of 
dwarf galaxy stars, if the internal mixing was not already 
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Figure 1. Compilation of stars with measured [C/Fe] and 
[Fe/H] in the stellar halo (squares), ultra-faint dwarf galax¬ 
ies (circles, hexagones, pentagones) and classical dSph galax¬ 
ies (triangles). [C/Fe] measurements are corrected to account 
for internal mixing processes (see text). CEMP-no stars are 
shown as filled symbols, upper limits with arrows. Stars 
with [C/Fe]> 0.7 and open symbols are CEMP-s/r stars. 
Filled gray squares are CEMP stars with no available mea¬ 
surements of r- or s-process elements. Colors/symbols iden¬ 
tify stars in dwarf galaxies with increasing total luminos¬ 
ity: fr om red to b l ue (s ee labels and text). References: Halo 


stars: IPlacco et al.l l| 20141) I Christlieb et al .l J2002); 


(2005)7 ~|Norris_et__al, d200jm [Caffai^t^al 


Frebel et al. 


(2014); [Hansen et al. (|2015l): I Bonifacio et al.N2015lh 
2015 b Segue 1: |N< 

Major II and Coma Berenice: 


Frebel et al 

et a lJ d201(t): IFrebel et al. I 11201 


Keller et al, 


Frebel et al.l ([2010 


Ursa 

Leo 


IV: [ Sim on et a ll d2010l b Bootes: lLai et a lJ (1201 lh; I Norris et alJ 
(l201ol): I Gilmo re et alj" d2013h . Dr a co: [ Cohen &; Hua nd (120 091); 


(2011). Carina: IVenn et al. 2012^). Sculptor: ( 

Frebel et al.l 2010aj; 

Tafelmever et al. 2010: Starkenbure: et al. 

201(3: Simon et al. 

2015; Kirbv et al. 20151; Jablonka et al. 20151; 1 Skuladottir et al. 


acc ounted for b y the authors, we used the online tool 
by IPlacco et alJ (2014) to self-consistently correct the 
[C/Fe] measurements. When not available from observa- 
t ions we si mply assumed [N/Fe] = 0.0, in agreement with 
IPlacco et al.l ( 2014 1. The data shown in Fig. 1 represents 
the largest sample of [C/Fe] measurements that have been 
homogeneously corrected to account for internal mixing 
processes. 

From Fig. 1 we can note that stars at [Fe/H]^ —4.5 
have only been found in the Galactic halo. Moreover, we 
can see that the frequency of CEMP-no stars among these 
“hyper-iron-poor stars” is extremely high: 8 out of 9 stars 
at [Fe/H]^ —4.5 are CEMP-no stars. These objects have 
peculiar chemical abundance patterns consistent with the 


yields predicted for primordial faint SN (e.g. llwamoto et alJ 
120051) . However, we also note that these extreme stars are 
part of a common, global trend, which involves stars in both 
the Galactic halo and dwarf galaxies. We see that on average 
[C/Fe] is higher in stars with lower [Fe/H], and it declines 
with [Fe/H]. The same trend affects the incidence of CEMP- 
no stars with respect to the overall stellar population. At 
[Fe/H]^ —3.0, the fraction of CEMP-no stars in the Galactic 
halo, (« 43%, Placco et al. 2014), is consistent with the 
overall fraction of CEMP-no stars in dwarf galaxies, where 
24 stars at [Fe/H]^ —3 have been observed, and 10 of them 
are found to be CEMP-no stars. This gives F C ehp(< —3) « 
42%. However, when we consider individual dwarf galaxies, 
the fractions are highly variable. 

The carbon measurements in the least luminous ultra¬ 
faint dwarf galaxies, log{L / Lq) < 4.0 (shown in Fig. 1 with 
red circles), are from high-resolution spectroscopic studies 
(see caption for references). Although less than 20 stars have 
been observed in these 4 systems, the overall frequency of 
CEMP-no stars is very high, and at [Fe/H]^ —3, 5 out 
of 6 stars are CEMP-no stars, Fcemp(^ —3) « 83%. With 
the only exception of seven star s in Bootes, one of which 
is a CEMP-no at [Fe/H]= —3.5 dGilmore et al.l [2013h . only 
low resolution spectroscopic studies are available for dwarf 
galaxies with luminosities between ultra-faint and classical 
dSph galaxies, e.g. Bootes, log{L / Lq) = 4.5±0.1, and Draco, 
log(L / Lq) = 5.45±0.10. Thus, slow and rapid n-capture el¬ 
ements have not been measured in these stars, preventing 
us from distinguishing between CEMP-no and CEMP-r/s 
stars. As discussed in Sec. 1, however, there is strong obser¬ 
vational evidence that CEMP-no stars should b e the d om- 
inant CEMP population at [Fe/H]< —3 (e.g. I Norris et alJ 
l2013ll . In these dwarf galaxies the fraction of CEMP (-no) 
stars at [Fe/H]^ —3 is high: 4 out of 6 stars in Bootes 
(« 66%), and 2 out of 5 stars in Draco 40%). The clas¬ 
sical dSph galaxies Carina and Sextans, log(L/ Lq) « 5.6, 
have been followed up at high-resolution, but only stars at 
[Fe/H]^ —3 have carbon measurements (see Fig. 1). 

In the Sculptor dSph galaxy, log{L / Lq) = 6.34 ± 0.16, 
many carbon me a surem ents are available from both 
low- (iKirbv et al.l 2015lj. an d high- re solution spectro¬ 
scopi c studies (IFreb el et al. _ 2010al: Tafelmever et_al] 

2010l: I starken burg et al~l [2013 : ISimon et al.l 120151 : 


Skuladottir et alJ 1201ol h In total, 10 stars have been 


found at [Fe/H]^ —3. After correcting for int ernal mixing, 
we fin d that the star Scl_ll_l_4296 observed by Simon et al. 
(l2015ll at [Fe/H]= —3.77 can be possibly identified as 
a CEMP-no star, [C/Fe]= 0.77 ± 0.34. This was not 
claimed by the authors, w ho used the luminosity depen¬ 
dent criteria developed by lAoki et alJ J2007t) to identify 
CEMP starfl .However, the carbon-excess is small and the 
errors big, making this classification very uncertain. Thus, 
the only reliable CEMP-no star discovered in Sculptor, 
[C/Fe]= 1.01 ± 0.1, has an unusua lly high iron-abundance, 
[Fe/H]« —2 ([Skuladottir et aT]l2015ii . which makes this star 
stand out with respect to the general trend in dwarf galaxies 
(Fig. 1). No CEMP-no stars have been found in Fornax, 
log(L/ Lq) ~ 7.25 ± 0.11, although carbon measurements 


1 There are small differences betw een the lAo ki et al.1 (120071 1 cri¬ 
teria and the carbon correction by lPlacco et al.l l20 1 4l . 
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are only available for stars at [Fe/H]> —2. For other 
ultra-faint dwarfs (e.g. Segue 2, Willmann, Hercules) or 
classical dSph galaxies (e.g. LeoT, Leol, LeoII, Ursa Minor) 
carbon measurements are not yet available. 

In conclusion, while Fcemp(^~ —3) in the Galactic halo 
is consistent with the overall fraction of CEMP-no stars at 
[Fe/H]^ —3 in dwarf galaxies, the individual Milky Way 
companions show that Fcemp(^ —3) strongly decreases when 
the luminosity of the galaxies increases. Does the fraction 
of CEMP-no stars depend on galaxy luminosity? And why 
does the only CEMP-no star observed in a more luminous 
dSph galaxy have an unusually high [Fe/H]« —2? Are these 
observational findings reconcilable with the idea that the 
birth environment of CEMP-no stars was polluted by pri¬ 
mordial faint SN, and thus consistent with what we see in 
the Galactic halo? 


3 COSMOLOGICAL MERGER-TREE MODEL 


We use the data-constrained cosmological code GAMETE 
(GAlaxy MErger Tree and Evolution, Salvadori^et_alJ 
2007ySalyadon^^emrral l200fll . 120121 : [Salvadori et al. 2014 


de Bennassuti et al. l2014h to link the properties of ancient 


stars in the Local Group with the early star formation and 
metal-enrichment processes. GAMETE describes the forma¬ 
tion of the Galaxy and its dwarf satellites in a ACDM frame¬ 
work, self-consistently accounting for the key physical pro¬ 
cesses driving the formation and evolution of high-redshift 
dwarf galaxies: (i) the transition from massive , and hence 
short-lived Population III (Pop III) stars, to normal Pop¬ 
ulation II (Pop II) stars; (ii) the gradual quenching of star 
formation in dwarf galaxies with increasing total masses due 
to the enhanced photo-dissociating and photo-ionizing radi¬ 
ation; (iii) the steady metal-enrichment of the Milky Way 
environment, or Galactic Medium , due to supernovae-driven 
outflows from star forming haloes. 

Our model is a statistical tool that enables us to 
study the most likely assembly and metal-enrichment 
histories of Local Group galaxies. The star formation and 
chemical evolution of present-day galaxies is traced across 
cosmic time by exploiting a sample of merging histories 
of a Milky Way-size dark matter halo ( Mmw = 10 12 Mq), 
reconstructed via a Monte Carlo algorithm based on the 
Exte n ded Press-Schechter form alism (see I Salvadori et al.l 
120071 : Ide Bennassuti et al.l I20l1 . for more details). Newly 
collapsed haloes are assumed to have a gas-to-dark mat¬ 
ter mass ratio equal to the baryonic cosmic fraction, 
Mg/M = Sib/£1 m and a chemical composition equal to 
the Milky Way environment at their formation epoch. 
Hence they are all primordial, Z = 0, before the onset of 
supernovae explosions. The main and new features of the 
model are summarized below, including the underlying key 
assumptions relevant for this work. We refer the reader to 
previous papers for more details. 


Star formation (SF). The SF is traced along the 
merger trees by adopting physically motivated hypotheses. 

• There exists a minimum halo mass to form stars, 
M s f(z), whose evolution accounts for the suppression of 
SF in progressively more massive objects due to the in¬ 
creasing photo-dissociating and photo-ionizing radiation 


(ISalvadori fe Ferrarall2009l . l2012l l. W hen the Milky Way en - 
vironment is fully reionized, z re i » 6 dSalvadori et alj|20l4 l. 
we assume that gas accretion is suppressed in haloes with 
virial temperatures T v i r ^ 2 x 10 4 K. 

• The SF rate, ip = e,M g /tff, which is regulated by 
the SF efficiency, e*, depends on the free-fall time, tff(z ), 
and mass of cold gas in each galaxy, whose gradual ac- 
cretion is described b y a numerically calibrated infall rate 
(ISalvadori et al .1120081 ). 

• In minihaloes with T v , r ;$ 10 4 K, the SF efficiency is as¬ 
sumed to be reduced as ch 2 = 2e«[l+ (T v i r / 2 x 10 4 if) ] _1 
to a ccount for the ineffective c ooling by molecular hydrogen, 
H 2 (ISalvadori fe Ferrarall2012l '). 

• Pop II stars with masses m = [0.1 — 100]Mq form 
according to a Larson Initial Mass Function, 'I'(to) = 

m~ 2.35 e -0.3SM e /m , jf the 

gas metallicity exceeds the crit¬ 
ical value, Z CI , which sets the minimal conditions for the 
formation of the first low-mass stars. This value can be ei¬ 
ther Z cr = (10 -4 — 10~ 3 ) Zq, if gas fragmenta t ion is driven 
by m etal-line cooling (e.g. iBromm et al.ll200ll : iFrebel et al.l 
M), or Z cr = (10 6 —10 4 )Zq, if it is due to thermal emis¬ 
sionbycoffisionally excited dust grains (e.g. fSchncidc r et al.l 
120021 : Omukai et a.1.1 f2005h . Following the most recent find¬ 
ings we set Z rT = 10~ 415 Z n (e.g. Caffau et al.l l201ll : 


ISchneider et al.ll2012l : ld~e Bennassuti et al . 2014l l. and we ex¬ 
plore the case Z CI = 10 Zq in Sec. 6. 

• Massive Pop III stars form if Z < Z cr . To work under 
the hypothesis t hat the early metal-enrichment is dominated 
by faint SN (e.g. [Cooke fe Madaull20l4 Ide Bennassuti et al.l 
12014 . we adopt the simplest prescription and assume that 
all Pop III stars have a characteristic stellar mass of 25 Mq 
and evolve as faint supernovae , which ex perience mixing and 
fallback fe.g. ISalvadori fc Ferrarall2012l ). 

Chemical enrichment. The contribution to the chem¬ 
ical enrichment by different stellar populations is traced 
in the code by accounting for the mass and metallicity 
dependent stellar lifetimes, a nd by including all chemi¬ 
cal elements from C to Zn (ISalvadori fe Ferraral l2012l l. 
For Pop III st ars evo l ving as faint supernovae we adopt 
the yie lds bv llwam oto et~ a l.l (l2005| l for the casfl Fol¬ 
lowing IWooslev fe Weaved ( 1995^1 we assume that stars 
with m > 40Mg do not contribute to metal enrichment, 
and we adopt their yields systematically halved in iron 
dTimmes et al.l Il995l l for massive 8 Mq ^ to ^ 40Mg 
stars that evolve as Type II supernovae (SNII), {Eshu ) = 
1.2 x 10 51 erg. For low and intermediate mass Asymptotic 
Gi ant Branch (AGB) stars, m < 8 Mp>, w e adopt the yields 
by Ivan den Hoek fe Groeneweeenl (Il997ll f or metallicities 
Z > 10~ 3 and by Mevnet fe Maeder ( 20021 1 for Z ^ 10~ 5 . 
All relevant equations describing the chemical enrichm ent o f 
star-forming galaxies can be found in lSalvadori et akl (l2008h . 

The contribution of supernovae Type la (SNIa) has 
been included by adopting the yiel ds a nd explosion energy 
(Esnia = 1.3 x 10 51 erg) bv llwamoto et all (|l999h , and the 
bimodal delay time distribution observationally derived by 


2 The total amount of Fe and C per stellar mass formed is 
equivalent to the integrated contribution of faint SN with m = 
(10 — 40 )M(T), as their yields are res caled with stellar masses (e.g. 
Fig.2 bv lde Bennassuti et al.1 (120141) ) 25 Mq and kinetic explosion 
energy £ P o P iii =0.7 x 10 51 erg 
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Figure 2. Predicted (contours) and observed (points) Fe-L re¬ 
lation for Milky Way dwarf galaxies. Contours identify regions 
containing the (68,95,99)% of the dSph candidates in 50 hier¬ 
archic al merger histories. Observed points are from I Kirby et alj 
teonh . Stars and labels underline individual dwarf galaxies. 




Mann ucci et al.ll (|2006h . At each time-step and for each star¬ 
forming halo of th e merger tree we c ompute the rate of SNIa 
by following iMatteucci et al.l (200d l and we fix the normal¬ 
ization constant, Asm a , to reproduce th e actua l rate of SNI a 
in the Milky Way, « (0.3/100)yr _1 (Canpellaro et alll999h . 

The chemical evolution of the gas is simultaneously 
traced in the star-forming haloes and in the surrounding 
Galactic Medium, or Milky Way environment, by includ¬ 
ing the effect of SN-driven outflows, which are regulated by 
the SN wind efficiency, e w . During SN explosions, metals 
present in the interstellar medium (ISM) are ejected out of 
the galaxy at a rate M’ff — ZMjjf where Z = Mz/M g is 
the metallicity of the gas, and M g 3 the gas ejection rate, 
which depends on the cumulative SN explosi on energy, a nd 
on the binding energy of the ho st halo (e.g. ISalvadori et al.l 
l200Sl: Ide Bennassuti et al]l20l4l . Heavy elements dispersed 
into the Milky Way environment are assumed to be instanta¬ 
neously mixed in this medium, so its total metallicity, Zgm, 
steadily increases across cosmic times. 

Universality of the free parameters. The free pa¬ 
rameters of the model are the star formation and SN-wind 
efficiency (e», e w ), the critical metallicity (Z cr ), and the 
fraction of stars that can give rise to SN la (Asnia). These 
are fixed to simultaneously reproduce the global properties 
of the Milky Wa y at z = 0, along with the MDF of Galactic 
halo stars (e.g. ISalvadori et al.l120071 : Ide Bennassuti et al.l 
2014). These free parameters are assumed to be the same 
for all the haloes of the merger tree and in all hierarchical 
merger histories of the Milky Way. Hence, the properties 
of different dwarf galaxies we will present in the following, 
are not the result of a fine tuning of the free parame¬ 
ters, but a consequence of the cosmological context in which 
we study the formation and evolution of these small systems. 

Selection and properties of satellite candidates. 

Galaxies that can survive as satellites of the Milky Way 
are selected among star-forming haloes of the merger 


trees that at any given redshift have dark matter masses 
that correspond to low-sigma density fluctuations, M < 
M 2 . 5 a(z). This assump tion is su pported by N-body simu¬ 
lations fe.g. iDiemand et al.l 120051 1. After selection, its sub- 
sequent evolution is follow ed in isolation down to 2 = 0 
(ISalvadori fc Ferrarall2009ll . 

In Fig. 2 we show the predicted Fe-Luminosity rela¬ 
tion for all selected candidates in 50 assembly histories of 
the Milky Way. We assume M/L = 1 to convert the stellar 
mass into total stellar luminosity, L * = M, x ( L/M ), and 
compare results with observations. We can see that model 
results match very well the observational data, including 
the nearly flat ([Fe/H]) value that is observed in ultra-faint 
dwarf galaxies. I 11 Fig. 2 the colored boxes identify galax¬ 
ies in different luminosity ranges, using the same color-code 
as in Fig. [T] These same colors will be used through the en¬ 
tire paper to distinguish among dwarf galaxies with different 
luminosities. 

In our model, ultra-faint dwarf galaxies are predicted 
to be H 2 -cooling minihaloes that form at 2 > 7.5, before 
photo-dissociation has suppress ed gas-cooling and star for - 
mation in these small systems (Salvadori fe Ferraral 12012 1. 
In contrast, classical dSph galaxies are predicted to form 
at later times through the merging of smaller progenitors. 
On average we find that the more luminous is a galaxy, the 
more massive is predicted to be its host halo, and th e lower 
its final assembly redshift (ISalvadori fe Ferrarall2009l f. More 
luminous galaxies, therefore, keep accreting gas from the 
Galactic Medium while it is increasingly polluted with the 
heavy elements ejected by low-mass star-forming galaxies. 

We finally underline that in our model we do not ac¬ 
count for mass transfer from binary companions. Thus, we 
can only investigate the incidence of CEMP(-no) stars that 
formed in carbon-enhanced environments. 


4 RESULTS: THE GLOBAL PICTURE 

We start by discussing the predictions of our cosmological 
model for the mean properties of dwarf galaxies in differ¬ 
ent luminosity ranges, which have been selected among all 
candidates in 50 different Milky Way realizations (Fig. [2]). 
The main results are summarized in Fig. [3] where we 
show, for dwarf galaxies with increasing luminosities (from 
top to bottom), predictions for the average: i) fraction 
of CEMP stars Fcemp (left), ii) normalized MDFs (mid¬ 
dle), and iii) low-Fe tails of the MDFs (right). The frac¬ 
tion of CEMP stars with respect to the total, Fcemp = 
Nl- 1 Efei A^([Fe/H])i/jV.([Fe/H])i, is calculated by av¬ 
eraging among the total number of dwarf galaxies within a 
given luminosity range, Nl, where A r * EMP ([Fe/H])i is the frac¬ 
tion of CEMP stars in the i-th galaxy, with a given [Fe/H] 
value. Similarly, we computed the mean MDFs by averag¬ 
ing among all 7V L dwarf galaxies in a given L range. For 
reference, Segue 1, Coma Berenice, and Ursa Major II have 
luminosities consistent with galaxies in the top panel (red), 
Bootes in the third (yellow), and Sculptor in the second from 
bottom (blue). 

Several interesting features can be noted in Fig. [3] (i) 
Independent of galaxy luminosity, we find that Fcemp = 1 
at [Fe/H]^ —5, and it rapidly decreases towards higher 
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Figure 3. Average properties of stars at different [Fe/H] in Milky Way dwarf galaxies with increasing luminosity (from top to bottom): 
L < 10 4 L 0 (red), 10 4 Lq < L < 10 4J3 L 0 (orange), 1O 4 ' 5 L 0 < L < 10 5 Lq (yellow), 1O 5 L 0 < L < 10 5 ' 5 L 0 (green), 1O 5 ' 5 L 0 < L < 
1O 6 L 0 (cyan), 10® L 0 < L < 1O 7 L 0 (blue), 10' Lq < L < 1O 8 L 0 (bottom, violet). Left panels: average fraction of CEMP stars (solid 
line) with ±1<t dispersion (shaded area). Middle panels: average MDF normalized to the total number of stars (histogram). Right panels: 
MDFs for stars at [Fe/H]< —3, showing both the total number of stars (filled histograms) and the number of CEMP-no stars (shaded 
gray histograms). The percentage shown is that of stars at [Fe/H]< —3 with respect to the total. 


[Fe/H], with a steeper decline in more luminous dwarf galax¬ 
ies. (ii) The shape of the normalized MDF dramatically 
changes with galaxy luminosity. In the least massive ultra¬ 
faint dwarf galaxies it is flat and covers a broad [Fe/H] range. 
As we move towards bigger galaxies, the MDF becomes more 
peaked, the peak is gradually shifted towards higher [Fe/H] 
values, and the low Fe-tail turns into a smaller fraction of 
the total, (iii) Stars with [Fe/H]< —3 are predicted to be 
found in all dwarf galaxies, but their relative contribution 
to the overall stellar populations strongly decreases when 
the luminosity of the galaxy increases, as shown by labels 
in the panels, (iv) On average, the cumulative number of 
stars at [Fe/H]< —4 is roughly of the same order of magni¬ 
tude, ps (1 — 2) x lO 3 A/ 0 , in all dwarf galaxies with the only 
exception of the faintest companions. 

All these features are simply a result of the hierarchi¬ 
cal galaxy formation process and the gradual metal enrich¬ 


ment and reionization of the environment out of which these 
galaxies form, which imprint their physical properties. In 
our scenario, the faintest dwarf galaxies are associated with 
H 2 -cooling minihaloes, which virialize from the Milky Way 
environment at 7.5 < z < 12, when —5 < [Fe/H]GM ^ 
—3. The chemical enrichment proceeds very smoothly in 
these small systems, which tra nsform gas into stars ver; 
inefficiently (Sec. 3 , see also Salv|rfoid&Ferrara 200! 
Webster et al.1 20141: Vincenzo et al.l 120141: I Romano et all 


I 


20151 : Bland-Hawth om et al.l 20151 ). Roughly a constant 


number of stars are formed at different evolutionary phases 
(or [Fe/H]), and the gas can reach high [Fe/H] before being 
either completely evacuated by the cumulative effect of SNe 
explosions, or photo-heated by the increa sing external ion¬ 
izing radiation (ISalvadori fc Ferrarall2012l ). This determines 
the characteristic shape of their MDF that is predicted, and 
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Figure 4. Predicted star formation rates of Milky Way dwarf 
galaxies in different luminosity ranges (colors as described in Fig. 
as a function of stellar ages (left) and redshift (right). 


observed, to be flat and to extend over a broad range of 
[Fe/H]. 

More luminous galaxies from via the assembly of these 
basics building blocks and more massive progenitors, which 
form afterwards. In this bottom-up scenario all dwarf galax¬ 
ies are expected to have, on average, some ultra-faint dwarfs 
among their parent haloes. Hence, they are predicted to 
share similar number of stars and MDF tails at the lowest 
[Fe/H]. However, the more luminous a galaxy, the less promi¬ 
nent is the contribution of the few thousand stars formed in 
these lowest-mass minihaloes into the final MDF. This is 
emphasized by the numbers reported in the right panels of 
Fig. [3] We can see that extremely metal-poor stars represent 
^ 50% of the stellar population in ultra-faint dwarfs, while 
they are only « 0.3% in dSph galaxies that have luminos¬ 
ity similar to Sculptor or higher. For this reason, the low-Fe 
tails seem to gradually “disappear” in the normalized MDFs 
as the luminosity of the galaxies increases. 

Since more massive dwarf galaxies complete their as¬ 
sembly at lower redshifts, they can continue accreting gas 
from a Milky Way environment that is increasingly pol¬ 
luted by heavy elements. The average [Fe/H]GM value at 
the formation epoch of the main halo, where the bulk of 
the stars form, sets the lower limit of [Fe/H] in the normal¬ 
ized MDF. Thus, on average, the more luminous a galaxy, 
the more pronounced is the shift of its MDF towards higher 
[Fe/H]. Furthermore, galaxies hosted by more massive dark 
matter haloes have bigger gas reservoir, and they can more 
efficiently convert gas into stars because of their higher 
T v j r > 10 4 K, which make them Lya cooling haloes (e.g. 
iMaio et al.ll2007h . This causes a gradually more pronounced 
peak in their MDFs. 


The fraction of CEMP stars in dwarf galaxies with dif¬ 
ferent luminosities is connected to all these previously pre¬ 
sented effects. Looking at the right panels of Fig. 02 we can 
see that IVcemp follows IV, at [Fe/H]< —4.5, in all dwarf galax¬ 
ies. These CEMP-no stars are predicted to form in (progen¬ 
itor) H 2 -cooling minihaloes that have been predominantly 
polluted by primordial faint SNe. Because of the small 
amount of iron (Vp e ~ 4 x 10 —7 ), and the huge amount of 
carbon (Y c fts 1CP 2 ) produced by faint SNe, long-lived Pop II 
stars can already start to form in these small systems[f] when 
the ISM is self-enriched up to [Fe/H] cr « —8, which corre¬ 
sponds to [C/Fe]« +4, and Z « Z cr = 10 ~ 4 Z@. After the 
formation of Pop II stars, normal SNII rapidly start to con¬ 
tribute to the chemical enrichment of both the ISM, gradu¬ 
ally decreasing the [C/Fe] level while [Fe/H] rises, and the 
Milky Way environment, increasing Zgm up to Z cr . In anal ¬ 
ogy to what is foun d in the Galactic halo (ICaffau et alJl201ll : 
iPlacco et aLll2014h we predict that carbon-normal stars can 
start to form when [Fe/H] U, — 4.7, which corresponds to 
Zgm > Z cr , and therefore to th e disappearance of Pop III 
stars (|de Bennassuti et al.ll20l3 ). 

Thus, we predict that Fcemp = 1 at [Fe/H] — 4.7 in 
all dwarf galaxies, and that this fraction rapidly decreases 
with increasing [Fe/H], because of the larger contribution 
of SNII in both self-enriched and newly formed galaxies. 
CEMP-no stars populating the MDF at —4 < [Fe/H]< —2, 
predominantly form in environment polluted by primordial 
faint SN and SNII. Thus they are all predicted to form in 
the very early stages of galaxy evolution. Interestingly, we 
can see that these CEMP-no stars efficiently form in dwarf 
galaxies like Bootes (yellow) and Draco (green), which are in 
the “transition region” between ultra-faint dwarfs and clas¬ 
sical dSph galaxies (Fig. [2j. This is because these galaxies 
are associated to the most massive among minihaloes, with 
T v ir « 10 4 K, which host Pop III stars as their least lumi¬ 
nous companion, but have higher star formation efficiencies 
(Sec. 3). More massive, classical dSph galaxies, have higher 
probabilities to have among their progenitor haloes Pop II 
galaxies that virialized from an enviro nment pre-enriched up 
to [Fe/HlGM A —4 by no rmal SNII l|Salvadori et, al.l[20071 : 
Ide Bennassuti et al.l 2014'). Carbon-normal stars efficiently 
form in these galaxies, causing the steeper decline of Fcemp 
with [Fe/H], 

In stars with higher iron-abundance, [Fe/H]> —2, we 
find that the carbon-enhancement might also come from 
low-metallicity (Z, ^ 10 ~ 3 Zq) AGB stars, which evolve on 
longer time-scales. This causes the small increase of the Fcemp 
values at [Fe/H]ss —1.5, which is visible in the left panels 
of Fig. 02 Note that in our model we do not account for bi¬ 
nary systems, so the carbon in these stars is not accreted 
from a companion, but it reflects the chemical composition 
of the birth environment. The products of AGB stars can 
be efficiently retained in the ISM of low mass galaxies after 
SNII have contributed to the chemical enrichment, and the 
galaxy can quietly evolve for a short period of time (e.g. 
ISalvadori fc Ferrarall2012lh The production of s-process el- 
ements from Z T 10“ 3 Zq AGB stars is still unclear (e.g. 
iFishlock et al.l [2014 ) and we do not account for it in our 


3 Where the “critical” iron-abundance is settled by the yields of 
faint SNe: [Fe/H] cr = log(Z CI /Zo) — log(Yz/Y Fe )+log(M 2 /M Fe )Q. 
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work. However, these stars are most likely expected to be 
CEMP-s stars, although their properties are probably dif¬ 
ferent from those formed in binary systems, which directly 
accrete material from an AGB companion. 

Finally, we should emphasize that, although Fcemp is 
steeper in more massive dwarfs, at [Fe/H]< —3 the results 
are all consistent within ±lcr errors. As already discussed, 
this is because these different dwarf galaxies are expected to 
share similar ancestors at high-redshift, namely low mass 
H 2 -cooling minihaloes , which represent the birth environ¬ 
ments of stars at [Fe/H]< —3. 

4.1 Star formation histories 

Another way to consider these findings is to look at the 
predicted average star formation histories (SFHs) of dwarf 
galaxies, which are shown in Fig. [4] as a function of stellar 
ages, and formation redshift. We can see that all dwarf galax¬ 
ies are predicted to host stars > 10 Gyr o ld, whic h is w hat 
is observed in the Local Group (e.g. ITolstov et al.ll2009n . At 
very high redshifts, z ^ 10 , we predict dwarf galaxies to 
share similar star formation rates, ip ~ 10 -3 — 10~ 4 Mq/ yr, 
independent of their total luminosity. It is during these cos¬ 
mic epochs that the low [Fe/H] tails of the MDFs are built 
up in the star-forming progenitor minihaloes. 

Ultra-faint dwarf galaxies are predicted to stop forming 
stars prior the end of reionization, z > 6 , because of heating 
by the external UV background that prevents gas cooling, 
and “sterilizes” those syst ems that still have some l eftover 
gas to fuel star formation dSalvadori fe Ferrarall201~2h . This 
is co nsistent w ith observations of SFH in ultra-faint dwarfs 
(e.g. lBrown et al in. Note that some of these “sterilized” 
minihaloes are predicted to evolve in isolation and survive 
until z « 2 , when they c an be observed as very-me tal-poor 
Damped Lya systems (ISalvadori fc Ferraral 120121 ). These 
systems might re-start forming stars at later times (e.g. 
iRicott i 2009; F aerman et alJl2013h and explain observations 
of Hi-rich dwarf gala xies that show recent episodes of star 
f ormati on ( e.g. Le oT: llrwin et al. I l2007l : Ide Jong et akllioosl : 
IWeisz et al.ll 20 li ~). 

More luminous galaxies continue to efficiently form 
stars to much later times, and hence have more extended 
and complex SFHs. In the right panels of Fig. [4]we see that 
ip is expected to gradually increase towards lower redshifts, 
reaching the maximum at the main assembling epoch of the 
host halo, and decreasing afterwards, when the galaxy is 
assumed to evolve in isolation. The SFH that has b een mea¬ 
sured in the Sculptor dSph galaxy dde Boer et al ] 120121 ) , is 
qualitatively consistent with the one that we predict for 
dwarf galaxies with similar luminosity (blue): it has a peak 
at ~ 13 Gyr, and then declines with cosmic time, lacking 
stars < 6 Gyr old. 

In this global picture we can therefore explain the vari¬ 
ation, with galaxy luminosity, of both the MDFs and star 
formation histories of dwarf galaxies. 


5 DATA COMPARISON 

We now test the results of our cosmological model against 
available observations for CEMP stars in nearby dwarf 
galaxies. We focus on three classes of dwarf galaxies that 


reside in different luminosity ranges, as shown in Fig. [2] 

i) the “classical” dSplr galaxy Sculptor, L = lO 6 , 3 ±o ' 2 Z/ 0 ; 

ii) the most luminous of the ultra-faint dwarf galaxies, 
Bootes, L = 10 4 5 ±o ' 2 Lq; and iii) the least luminous ultra¬ 
faint dwarfs, L < 10 4 Lq: the combination of Seguel, Coma 
Berenice, Ursa Major II, and LeoIV. As pointed out in Sec. 2, 
the available data for Bootes are mostly blind to CEMP sub¬ 
classes (s, r, and no). This represents a possible caveat for 
our comparison with CEMP-no models. Still, most of the 
CEMP stars found in this system have [Fe/H]< —3, which 
is t h e typ ical [Fe/H] range of CEMP-no stars (e.g. Fig. 1 of 
I Norris et al .11 20131 ). 

Given the low number of [C/Fe] measurements in dwarf 
galaxies (see Fig. [5^,, [ 6 ^,d), we com puted the uncertainty 
of Fcemp by using the results of ICehrelsl (jl986i ) for Poisson 
statistics. We derived the 1 <t upper (lower) confidence limits 
as ^cemp = 1.841/A* ([Fe/H]) in case of non-detection (F/emp = 
0.0), and as F^ = 3.300/A*([Fe/H]) for single detection 
(F/S£ = 0.173/A*([Fe/H])), where A*([Fe/H]) represents the 
total number of stars with available carbon measurements 
in different [Fe/H] ranges. So the fewer the measurements, 
the larger the uncertainties. 

We used a Monte Carlo technique to randomly select 
from the mean theoretical MDFs, a number of stars equal to 
the total number of [Fe/H] measurements in different dwarf 
galaxies, Ntot■ We then constructed the average ±lcr errors 
on the MDF by iterating this procedure 100 times. These 
errors, which are shown in Fig. [5]r and Fig. [ 6 ] 3 ,c together 
with the dispersion among different dwarf galaxies, provide 
an estimate of the data-driven uncertainties, i.e. they are 
the errors we always have to deal with while comparing the¬ 
oretical models with small numbers of measurements. 

5.1 Carbon-enhanced stars in Sculptor 

The predicted average properties of Sculptor-like dwarf 
galaxies are shown in Fig. [5] We note that there is only one 
det ection of a CEMP-no st ar in Sculptor at —2.5 <[Fe/H]^ 
—2 ([Skuladottir et al.t2015l) . from which we get F^ p = 0.06. 
In the other [Fe/H] bins, with zero detections, F^ p is simply 
set by the number of [C/Fe] measurements. So the decline of 
-Fcemp towards higher [Fe/H] is simply a consequence of the 
larger number of stars observed at increasing [Fe/H]. 

From Fig.[5]i we can see that the theoretical predictions 
are consistent with F/^Jp. As we already discussed in Sec. 4, 
because of the early chemical enrichment by primordial faint 
SN, Fcemp rapidly declines with [Fe/H], ranging from unity 
at [Fe/H]^ —4.75 down to values < 0.005, at [Fe/H]> —1. 
Carbon-normal stars begin to form when [Fe/H]^ —4.7, and 
in Sculptor they are predicted to be the majority of the 
stellar population (> 50%) already at [Fe/H]is —4.25. 

This result has two implications: i) to catch second- 
generation stars imprinted mainly by primordial faint SNe 
in Sculptor, we should look among stars at [Fe/H]< —4.75, 
as in the Galactic halo; ii) to increase the probability of find¬ 
ing CEMP stars in Sculptor we should follow up [Fe/H]< —3 
stars, for which we predict Fcemp ^ 0.065. Our ability to 
find these CEMP stars will remain naturally limited by the 
absolute number of stars that exist at these low [Fe/H]. 
Our model predicts, and this is supported by observations, 
that stars at [Fe/H]< —3 are intrinsically rare in luminous, 
Sculptor-like galaxies, only representing < 3% of the to- 
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Figure 5. Predictions vs observations for the Sculptor dSph 
galaxy. Top panel: the fraction of CEMP stars in different [Fe/H] 
bins. The solid line shows the average value in the model, and 
the shaded area the Tier dispersion among different Sculptor¬ 
like dwarf galaxies (as Fig. CD- Upside down triangles are upper 
limits for the observed fraction of CEMP stars based on avail¬ 
able data. Labels indicate the number of [C/Fe] measurements 
in each [Fe/H] bin (references in Fig. [TJ. Middle panel: the ob¬ 
served (points with errorbars) and predicted MDF (histograms 
with shaded uncertainty). We show the Tier dispersion among 
Sculptor-like dwarf galaxies in 50 Milky Way possible assembling 
histories (light gray), and among 100 Monte Carlo sampling of 
the average MDF to the number of stars observed (dark gray). 
Labels indicate the number measurement s in each [Fe/H] bin 
(iKirbv et alj|20lH : IStarkenburg et alj|2013h . Bottom panel: con¬ 
ditional probability to observe a star at a given [Fe/H] that it is 
also carbon-enhanced. Labels indicate the number of [C/Fe]> 0.7 
measurements in each [Fe/H]. 

tal stellar population (see Fig. 0. This point is highlighted 
in Fig. 0), where we compare the predicted and observed 
MDFs, that are normalized to the total number of stars. 
Noticeably, stars at [Fe/H]< —3 are almost invisible in the 
normalized function, which is dominated by more Fe-rich 
stars. Clearly, the paucity of stars at this low [Fe/H] makes 
it very challenging to search for CEMP-no stars in Sculptor, 
where only the brighter stars at the tip of the RGB can be 
followed up. Between —3.5 <[Fe/H]< —3, for example, we 
should roughly double the number of stars found in Sculptor 
to be able to observe one CEMP-no star among them. 

Thus, we can ask a different question: what is the joint 
probability to observe a star that has a given [Fe/H] value, 
and that is also carbon-enhanced? Fig. 0 shows this joint 
probability distribution. This has been computed by com¬ 
bining the two independent functions: the fraction of CEMP 
stars predicted by the model, and the observed MDF normal¬ 
ized to the total number of stars, P 0 bs = Fcemp x N*/Ntot- 
The probability, P 0 b 3 , is maximal at [Fe/H]~ —2, where we 
predict the highest absolute number of CEMP stars. This 
result naturally explains why the first carbon-enhanced star 
in the Sculptor dSphs has been serendipitously found at such 
a high [Fe/H], In conclusion, the observed MDF, the de¬ 
rived upper limits for Fcemp, and the iron-abundance of the 
only CEMP star detected in Sculptor, are all consistent with 
the hypothesis that faint SNe were the main contributors to 
primordial metal enrichment. 


5.2 Carbon-enhanced stars in ultra-faint dwarfs 

We now explore the predictions of our model for less lu¬ 
minous ultra-faint dwarf galaxies, where many CEMP stars 
have been found (Fig. 0. In Fig. 0 we show the results for 
Bootes-like dwarf galaxies and for the least luminous ultra¬ 
faint dwarfs with L < 10 a Lq. As for Sculptor we select these 
galaxies on the bases of their total luminosity among all the 
available candidates predicted by the model (Fig. [2]). 

By inspecting Fig. [6jr and Fig. 01, we can see that be¬ 
cause of the low total number of stars observed in these small 
systems (« 40 in Bootes, and fts 15 in the faintest ultra-faint 
dwarf galaxies combined all together) the upper/lower lim¬ 
its on Fcemps are extremely high/low, and hence naturally 
consistent with model predictions. In these ancient systems, 
the fraction of RGB stars with respect to the total is esti- 
matecQ to be Argb/A* 0 * « 0.001. This implies that in both 
cases > 50% of RGB stars have been already followed-ujjf| 

We note that in the model, the dispersion among differ¬ 
ent galaxies is huge in the case of the faintest dwarf galaxies 
(Fig. 01), and at [Fe/H]= —3.5 the Tier error is consis¬ 
tent with Fcemp « 0.8. This is because these faint galaxies 
are associated with the least massive H 2 -cooling minihaloes 
that are able to form stars, and which evolve in isolation. 
These systems can either virialize from a primordial birth- 
environment, Z < Z cr , and hence Fcemp ~ 1 for [Fe/H]< —3 
because they host primordial faint SN, or from a medium 
that has been pre-enriched up to Z > Z cr by the prod¬ 
ucts of SNII. In the latter case they will only host Pop II 
stars, and hence Fcemp ~ 0. These different formation paths 
cause the large spread in Fcemp at [Fe/H]< —3. According 
to our model, the large fraction of CEMP-no stars observed 
in the faintest dwarfs at [Fe/H]< —3 (« 80%) suggests that 
these systems are truly Pop III galaxies, which have experi¬ 
enced primordial star formation. However, we clearly need 
better measurement statistics. A quest for more data also 
emerges while comparing the predicted and observed MDFs 
(Fig. 0:). In the faintest dwarf galaxies, the errors induced 
by the low number of measurements (15 stars) are larger 
than the dispersion of the model among different dwarf 
galaxies and Milky Way merger histories. In contrast, these 
errors are equal (lower) than the intrinsic model uncertain¬ 
ties in Bootes (Sculptor), where ~ 40 stars (« 700) have 
been observed. 

Fig-0? f show the joint probability to observe a star in a 
given [Fe/H] range, that is also carbon-enhanced. If we com¬ 
pare this function with the one derived for Sculptor (Fig.0), 
we note a clear trend: the fainter the galaxy is, the higher is 
the overall probability to observe a star that is also carbon- 
enhanced, which is the integral of P 0 t, s over the [Fe/H] range. 
This explains why it is much easier to find carbon-enhanced 
stars in these small systems. Moreover, P 0 b s is maximum 
at lower [Fe/H] in less luminous dwarf galaxies. For Bootes 
we End that P 0 b s is maximum at —4 ^[Fe/H]^ —2.5, while 
for the faintest dwarfs at [Fe/H]sC —3. In both cases the pre- 

4 We used PARSEC isochrones llBressan et al 1 l2012l : IChen et al. 
|2014| ; iTang et alj l2014lh and the CMD generator available at 
http://stev.oapd.inaf.it/cmd 

5 For the ultra-faint dwarfs this only is true if we consider the 
combination of the four systems at L < 10 ^Lq. In Segue I all 
RGB stars have already been followed up. 
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Figure 6. Same as Fig.[5]but for Bootes-like dwarf galaxies (left), and for the least luminous ultra-faint dwarfs with L < KpLg (right). 
Measurements are the same reported in Fig. [T] For Fcemp we show both upper and lower limits from observations (see text). 


dieted [Fe/H] ranges coincide with [Fe/H] values of observed 
CEMP-no stars in these small systems. 

Finally, we should note that for the faintest dwarf galax¬ 
ies the model predicts a non-negligible fraction of CEMP 
stars at [Fe/H]> —2 (upper panel). As explained in Sec. 4, 
these stars are predicted to form in environment polluted by 
the products of AGB stars with Z < 10~ 3 Zq, which can be 
retained by these small galaxies when SN have already ex¬ 
ploded. Interestingly, a CEMP-s star at [Fe/H]« —1.5 that 
does not show evidence for a bina ry companio n has been re¬ 
cently detected in Segue 1 (lFrebel et al.ll2014r l. So, although 
more data are required to solidly assess if this star is in a 
binary system, this observation supports the idea that the 
physical mechanism we propose may contribute to the for¬ 
mation of CEMP-s in these small systems. 


6 THE LOW-[FE/H] TAIL OF SCULPTOR 

The results of our cosmological model show that CEMP-no 
stars, similar to those found in the Galactic halo and ultra¬ 
faint dwarfs, are not necessarily missing in classical dSph 
galaxies, such as Sculptor, but they could be hidden among 
rare [Fe/H]< —3 stars, which represent < 3% of the overall 
stellar population. The possible existence of a CEMP-no star 
at [Fe/H]= —3.77 (see Section 2) might be the first confirma¬ 
tion that this is indeed the case. Sculptor-like dSph galaxies, 
furthermore, are predicted to be among the best systems to 
find the most iron-poor stars, because their MDFs might ex¬ 
tend down to [Fe/H]< —5, if faint SN dominated the early 
metal enrichment Fig. [3] The key question is: how much 
do we need to enlarge the current stellar sample to unveil 
the lowest-Fe tail of the Sculptor MDF and catch the most 
pristine CEMP stars? 

Fig. [3 shows how many stars at [Fe/H]< —3 are pre¬ 
dicted to emerge in Sculptor-like dwarf galaxies by increas¬ 
ing the total number of [Fe/H] measurements. Fig. ex¬ 
hibits results for the actual number of observed stars, « 700 
stars, which roughly correspond to « 25% of all RGB stars 
in Sculptor. Model predictions agree quite well with obser¬ 
vations. Fig. [3 j shows how the low-Fe tail of Sculptor will 


appear by following up all stars down to magnitudes V ^ 20, 
which corresponds to « 16 00 stars in the observed Color 
Magnitude Diagram dde Boer et al.lEoi2h . We find that in 
this case we might be able to observe 12 ± 8 stars with 
[Fe/H]< —4, the (40 ± 20)% of which are predicted to be 
CEMP-no stars. These observations can be achieved with 
present-day instrumentation (e.g. ESO VLT with X-Shooter, 
FLAMES, and UVES, or Keck/DIMOS) to add crucial in¬ 
formation to our understanding of the properties of first 
sta rs and early galax y formation processes. See for exam¬ 
ple (Kirb^etjil] ( 20091) . who used Keck/DIMOS to obtain 
medium-resolution spectra (6400 — 9000 A, R « 6500) for 
V ss 20 stars in the central region of Sculptor and derive 
their C and Fe abundances. Future generation of telescopes 
will allow us to measure [Fe/H] for stars below the main se¬ 
quence turn-off in nearby galaxies, like Sculptor. For exam¬ 
ple, the limit given for MOSAIC spectroscopy on the ESO- 
ELT (Evans et al. 2015) is V « 25 at R as 15, 000 - 20, 000. 
This will dramatically increase the number of stars that can 
be observed at high spectral resolution in nearby galax¬ 
ies. The total number of RGB and MS stars in Sculptor 
with V ^ 23 is sa 20000, and from these we can expect 
80 ± 22 stars at [Fe/H]^ —4 (Fig. [T}:) , and 16 ± 10 stars 
with [Fe/H]< —4.7, where the fraction of CEMP-no stars 
is expected to be 100%. With these observations, further¬ 
more, we might be able to accurately con strain t h e cri ti¬ 
cal m etallicity value. As already shown in ISalvadori et al.l 
(l2007ll (Fig. 7), decreasing Z cr in the most feasible range, 
Z cr ~ 10 -4 — 10 -6 Zq, mainly affects the number of long- 
lived stars at [Fe/H]< —4. Thus, the differences between 
Z cr models become observable in Sculptor when a signifi¬ 
cant number of stars are followed-up, and the lowest-Fe tail 
starts to emerge (Fig. 7c). 


7 SUMMARY AND DISCUSSION 

We used a robust, data-constrained merger-tree model for 
the possible formation histories of the Milky Way and 
its dwarf satellites to predict the fr equency of CEMP-no 
stars in nearby dwarf galaxies (e.g. I Salvadori et al.l 120081 : 
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Figure 7. Number of stars at [Fe/H]< —3 that can be observed in 
Sculptor by increasing the sample of [Fe/H] measurements. From 
top to bottom we show results for: i) the current statistics (~ 700 
RGB stars), where the available data are shown as points with 
Poissonian errorbars (Fig. 5 for references); ii) stars with V ^ 20; 
iii) stars with V ^ 23; iv) all stars. Shaded area show the dzlcr 
errors obtained from the Monte Carlo random selection technique 
(see text). Solid histograms with shaded area are results for our 
fiducial model, Z cr = 10 —4-15 Zq. Red dashed histograms show 
the same results for the case Z cr = 10 ~ 6 Zq. 


ISalvadori fc Ferraral l200fll . 120121 : ISalvadori et al.l 120141 '). We 
have shown that the model can successfully explain the vari¬ 
ation of the average MDFs and SFHs observed in dwarf 
galaxies with incr easing luminosities fe.g. lKirbv et al.ll201ll : 
IWeisz et al.l 120141 1 by accounting for star-formation in H 2 - 
cooling mini-haloes, M = 10 6 ’ 5 — 10 8 M©, with a self- 
consistent treatment of the reionization and metal enrich¬ 
ment of the Milky Way environment. 

By assuming that primordial faint SN, with mixing 
and fallback, dominated the early chemical enrichment, 
as suggested by observations of Galactic halo stars (e.g. 
Iwamoto_e^_aL|_ 2005|^ Cooke fc Madaul 120141 ; iMarassi et alJ 


20 111 : de Bennassuti et alT 2014 1. we analyze the average 


properties of dwarf galaxies in different luminosity ranges, 
and show that: 


• CEMP-no stars should exist in all dwarf galaxies within 
the observed luminosity range, 10 2 ' 5 L© < L < 10' 5 Z/©. 

• Independent of galaxy luminosity the relative fraction 
of CEMP-no stars increases towards lower [Fe/H], reaching 
100% at [Fe/H] < - 4.7. 

• As the galaxy luminosity increases, the overall proba¬ 
bility to observe CEMP-no stars decreases, and the [Fe/H] 
range in which they are most likely to be found is shifted 
towards higher values. 

• In classical Sculptor-like dSph galaxies, the probability 
to find CEMP-no stars is almost an order of magnitude lower 
than in the faintest Milky Way companions, L < 10 4 L©, and 
it is maximal, P 0 b s = 0.02 at [Fe/H]~ —2. 


The results explain why it is easier to discover CEMP- 
no stars in ultra-faint d warf galaxies than in more lumi¬ 
nous classical dSp hs ('e.g. iNorris et al1l20ld : lLai et al.ll201ll : 
iFrebel et al.ll2014h . and also why the only CEMP-no star ob¬ 
served i n Sculptor has been fou nd at an unexpectedly high 
[Fe/H] (ISkuladottir et alJ 120151') . In particular, our model 
shows that CEMP-no stars at [Fe/H] ^ — 4.7 form in en¬ 
vironments predominantly polluted by primordial faint SN, 
which therefore have extremely high [C/Fe] values. On the 
other hand, those at —4.7 [Fe/H] —2.0 are imprinted by 
both primordial faint SN and low metallicity, Z < 10 ~ 3 Zq, 
SNII, in agreement with recent findings by Bonifacio et al. 
(2015). At [Fe/H] — 2.0, we find that CEMP stars can 
also form in an ISM polluted by the products of AGB stars 
with Z < 10" - 3 z 0 . Hence, they may potentially be enriched 
by s-process neutron capture elements (i.e. CEMP-s), which 
are not included in our model. 

Our findings for CEMP-no stars, are a consequence of 
both the extremely low Fe-production, and high C, from 
faint SN, and of the cosmological context in which the hier¬ 
archical assembly of dwarf galaxies occurs. H 2 -cooling mini¬ 
haloes, or ultra-faint dwarf galaxies, L < 10 5 L©, are thus 
predicted to be the high-redshift progenitors of more mas¬ 
sive “classical” dSph galaxies, and the environment of for¬ 
mation for stars with [Fe/H]< —3. The low-Fe tails of dwarf 
galaxies are thus expected to form in these common build- 
ing blocks, some of wh ich experienced Pop III star formation 
(ISalvadori et alj|2014j l. Such low-Fe tails can extend down 
to [Fe/H]< —4.7 if built-up upon the chemical products of 
primordial faint SN. 

As the galaxy luminosity increases our model shows 
that the average MDFs become more peaked , and shifted 
towards higher [Fe/H] values, as is observed. Thus, stars 
at [Fe/H]< —3 become a lower fraction of the total stel¬ 
lar populations: from ^ 40% for ultra-faint dwarf galaxies, 
L < 10 s L©, to < 0.2% for L > 10'L©. This is why in more 
massive dwarf galaxies the overall probability to observe a 
CEMP-no star decreases, and the [Fe/H] range in which 
they are more likely to be found increases. In Sculptor, L ~ 
10 6,3 L©, we find that the fraction of CEMP-no stars mono- 
tonically decreases from Fcemp ~ 1 at [Fe/H]^ —4.75 down 
to Fcemp < 0.005 at [Fe/H]^ —1. However, the probability to 
observe one of these stars is higher at —2.5 <[Fe/H]^ —2 
because more stars can be observed in this range. 

A key prediction of our work is that galaxies with dif¬ 
ferent luminosity are expected to share, on average, similar 
MDF tails and fractions of CEMP-no stars at the lowest 
[Fe/H]. At the moment, this hypotheses is supported by the 
observations of a very similar fraction of CEMP stars at 
[Fe/H]< —3 in the Galactic halo and in the overall sample 
of nearby dwarf galaxies, Fcemp(^ —3) ~ (43 — 42)%. How¬ 
ever, more data need to be collected for the Milky Way com¬ 
panions. This also emerges from our analysis of the model 
uncertainties, which in the faintest dwarf galaxies are dom¬ 
inated by the low number statistics. 

Finally, we should be clear that our model cannot ex¬ 
clude alternative scenarios for the formation of CEMP-no 
stars, suc h as metal pollution by massive rotating primordial 
stars ('e.g. iMevnet et ahll2006l : iMaeder et al .1120151 ). Further¬ 
more, we did not explore the effects of different primordial 
Initial Mass Functions on the fraction of CEMP-no stars in 
dwarf galaxies. Instead, we simply worked under the hypoth- 
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esis that Pop III stars evolve as faint SN with mixing and 
fallback, since these pristine stars should dominate the early 
chemical enrichment to successfully explain observations o f 
CEMP-no stars in the Galactic halo (e.g. ICooke fe Madaul 
120141 ; Ide Bennassuti et al. inn. Our findings show that to 
further test the predominant role of primordial faint SN in 
the early Universe we need to at least double the current 
stellar sample in Sculptor, by measuring [Fe/H] and [C/Fe] 
in stars in different region of the galaxy down to magni¬ 
tudes V ^ 20. Deeper observations (V Sj 23) should enable 
the discovery of peculiar CEMP-no stars in Sculptor simi¬ 
lar to those found in the Galactic halo at [Fe/H]< —4. This 
will provide crucial information not only on the nature of 
first stars and critical metallicity value, but also on the for¬ 
mation of the Galactic halo, and hence on the underlying 
hierarchical galaxy formation process. 
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